Wigner function's negativity demystified 
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As early as 1932 Wigner denned the joint distribution for the coordinate and momentum of a 
quantum particle. Despite a drawback of being sometimes negative, the Wigner distribution has 
stood the test of time and found many applications. Having demonstrated that the Wigner function 
of a pure quantum state is a wave function in the specially tuned Dirac bra-ket formalism, we argue 
that the Wigner function is in fact a probability amplitude for the quantum particle to be at a 
certain point of the classical phase space. Since probability amplitude need not be positive, our 
findings elucidate the long-standing mystery of the Wigner function's negativity. Additionally, we 
establish that in the classical limit, the Wigner function transforms into a classical Koopman-von 
Neumann wave function rather than into a classical probability distribution. As a result, contrary 
to widespread beliefs, the volume of negative regions in the Wigner distribution cannot quantify the 
degree of quantum character. 



In his seminal work [1 j, Wigner denned the combined 
distribution of the quantum particle's coordinate and mo- 
mentum in terms of the wave function. Since then, this 
function, bearing his name, plays a paramount role in 
the phase space formulation of quantum mechanics [2] , is 
a standard tool for establishing the quantum-to-classical 
interface |3j H] , and has a broad range of applications in 
optics and signal processing [5]. Techniques for the ex- 
perimental measurement of the Wigner function are also 
developed [6HH]. Despite its ubiquity, the Wigner func- 
tion is haunted by the obscure feature of possibly being 
negative. Wigner [9] upgraded this peculiarity to a mys- 
tery by demonstrating that his function is the only one 
satisfying a reasonable set of axioms for a joint proba- 
bility distribution. Subsequently, the Wigner function's 
negativity motivated the development of a mathematical 
framework for handling negative probabilities [lOj [TT] . 
Moreover, it is widely believed that the presence of neg- 
ative regions in the Wigner distribution is an authentic 
sign of quantum character. Even though doubts were 
cast on the validity of such an interpretation [12 j, the 
significance of the negativity has remained unclear. 

In this Letter we clarify the origin of the negativity 
by advocating the following interpretation: The Wigner 
function is a probability amplitude for a quantum parti- 
cle to be at a certain point of the classical phase space, 
i.e., the Wigner function is a wave function analogous to 
the Koopman von-Neumann (KvN) wave function of a 
classical particle. 

Around the time the Wigner distribution was con- 
cieved, Koopman and von Neumann [13-16 recast classi- 
cal mechanics in a form similar to quantum mechanics by 
introducing classical complex valued wave functions and 
representing associated physical observables by means of 
commuting self-adjoint operators. In particular, it was 
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postulated that the wave function \^(t)} of a classical 
particle obeys the following equation of motion: 

<^ !*(*)> = £ !*(*)>> l = ^\ x -u\x)\ p , (i) 
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[x,p] = [x,X p ] = \p,X x ] = [XxAp] = 0. 



(2) 



Without loss of generality one-dimensional systems are 
considered throughout. Since the self- adjoint opera- 
tors representing the classical observables of coordi- 
nate x and momentum p commute, they share a com- 
mon set of orthogonal eigenvectors \px) such that 1 = 
J dpdx \px) (px\. In the KvN classical mechanics, all 
observables are functions of x and p. The expectation 
value of an observable F = F(x,p) at time t equals to 
(V(t)\F | The probability amplitude (px\^(t)) for 

a classical particle to be at point x with momentum p at 
time t is found to satisfy 
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(px\^(t)) = 0. 



(3) 



This is the evolution equation for the classical wave func- 
tion in the xp- representation, where x = x, X x = 
—id/dx, p = p, and X p = —id /dp in order to satisfy 
the commutation relations Utilizing the chain rule 
and equation ([3|, we obtain the well known classical Liou- 
ville equation for the phase-space probability distribution 
p(x,p;t) = \(px\V(t))\ 2 , 
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p(x,p;t) = 0. 



(4) 



Note that the classical wave function and the classical 
probability distribution satisfy the same dynamical equa- 
tion, which reflects the physical irrelevance of the phase 
of a classical wave function. 

In recent work [17 , we put forth the Ehrenfest quanti- 
zation as a systematic theoretical framework for deducing 
equations of motion from the evolution of average values. 
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First, starting from the Ehrenfest theorems, we obtained 
the Schrodinger equation if the momentum and coordi- 
nate operators obeyed the canonical commutation rela- 
tion, and the KvN equation if the momentum and co- 
ordinate operators commuted. Then, applying the same 
technique to the Ehrenfest theorems, 

mj t <¥„(t)| x q |¥„(t)) = <¥„(*) | p, |¥„(t)) , 

j t (V K (t)\p q |*„(*)> = <*„(*)! - £/'(£,) |* K (t)> , (5) 

with a generalization = zHk (0 ^ « ^ 1) and 

demanding a smooth classical limit n — » 0, we have shown 
that 

^ i*„(t)> = ?i qc i*„(t)> , 

^ c = m PK + -. U \ X - ~2 p )-- K U { X+ Y X o)> 

X q = X — HkX p /2, p q =p + hK,\ x /2, (6) 

where x q and p q represent the quantum coordinate and 
momentum respectively, x, p : \ x , and X p are the same 
classical operators as in equation ([2|, and k denotes the 
degree of quantumness/commutativity: n — » 1 corre- 
sponds to quantum mechanics whereas ft — >• recovers 
classical mechanics, lim^o Wqc = fit. See figure [l] for 
a pictorial summary of these derivations thoroughly pre- 
sented in [T7] . 

A crucial point for our current discussion is that this 
unified wave function \if? K ) (t is dropped henceforth) in 
the xp-representation is proportional to the Wigner func- 
tion W, 

(px \Sfr K ) = y/27rhKW(x,p), 
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where p K denotes the density matrix with h replaced by 
Hk [17]. Straightforward calculations show that the nor- 
malization condition for the unified wave function implies 
that the density matrix must correspond to a pure state, 

<* K |*0 = l<=>fc = p K . (8) 



Setting p K = \(j) K ) ((/) K \, we can demonstrate that all ex- 
pectation values for the state \<p K ) coincide with those for 

i*o, 



(* K | G(x q )F(p q ) |* K ) = (<p K \ G(x q )F(p q ) |0O , (9) 



where G(x q ) and F{p q ) are arbitrary functions of the 
quantum position and momentum, respectively. 

Equations ([8| and ([9| reveal that the Wigner distribu- 
tion of a pure state indeed behaves like a wave function. 
As shown in figure [I] the Wigner function's dynamical 
equation (|6| transforms into the evolution equation for a 
classical KvN wave function ([!]). Hence, in the classical 
limit, the Wigner function maps a quantum wave func- 
tion into a corresponding KvN classical wave function 
rather than a classical phase space distribution. Since the 
vectors \px) are identical in both KvN and Wigner rep- 
resentations, W(x,p) is proportional to the probability 
amplitude that a quantum particle is located at a point 
(x,p) of the classical phase space. Note it is important 
to distinguish the classical (x,p) and quantum (x q ,p q ) 
phase spaces because the notion of a phase-space point 
arises naturally only in the commutative classical vari- 
ables (x,p). One may take this distinction further and 
interpret the Wigner function as the KvN wave function 
of a classical counterpart of the analogous quantum sys- 
tem. In any case, the Wigner function, as any other wave 
function, need not be strictly positive, and its negativity 
does not necessarily imply the quantum nature of either 
the system or dynamics. 

The following statement, proven in [17], serves as an 
illustration of the interpretation above: The equality 
H qc = HL holds solely for a quadratic potential; namely, 
quantum and classical harmonic oscillators with the iden- 
tical spring constant share the same eigenstates. Hence, 
the Wigner and KvN wave functions of these eigenstates 
coincide. 
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FIG. 1: The derivation of quantum mechanics, classical mechanics, and quantum mechanical phase space representation within 
the Ehrenfest quantization proposed in [17] . 
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